:arXiv:chao-dyn/9605015vl " 31 May 1996; 



Anomalous scaling in random shell models for 

passive scalars 

A. Wirth 1 & L. Biferale 2 
February 5, 2008 

1 CNRS, Observatoire de Nice, B.P. 229, 06304 Nice Cedex 4, France. 
2 Universita di Roma "Tor Vergata" , Dip. di Fisica 
Via della Ricerca Scientifica 1, 00133 Roma, Italy. 



Abstract 

A shell-model version of Kraichnan's (1994 Phys. Rev. Lett. 72, 1016) passive scalar 
problem is introduced which is inspired from the model of Jensen, Paladin and Vulpiani 
(1992 Phys. Rev. A45, 7214). As in the original problem, the prescribed random velocity 
field is Gaussian, delta-correlated in time and has a power-law spectrum oc k^, where k m is 
the wavenumber. Deterministic differential equations for second and fourth-order moments 
are obtained and then solved numerically. The second-order structure function of the passive 
scalar has normal scaling, while the fourth-order structure function has anomalous scaling. 
For £ = 2/3 the anomalous scaling exponents Cp are determined for structure functions up to 
p = 16 by Monte Carlo simulations of the random shell model, using a stochastic differential 
equation scheme, validated by comparison with the results obtained for the second and fourth- 
order structure functions. 

ROM2F/96/29 



1 Introduction 

A most striking property of a passive scalar quantity advected by a fully developed turbulent flow 
is its spatio-temporal intermittent behavior (see Ref. [|IJ for a recent review). The time evolution 
of the scalar field 9 is described by the partial differential equation : 

d t e + (v • v)0 = kA6 + /, (i) 



where v is the advecting velocity field, k is the molecular diffusivity and / is an external forcing. 



Intermittency in passive scalars advected by highly chaotic flows is thought to be connected to 
the existence of self-similar processes transferring fluctuations from large to small scales. 

Self-similarity is observed in the power-law behavior of structure functions, S p (r), in the inertial 
range, that is moments of ^-increments at scales where neither external forcing nor molecular 
damping are acting : 

S p {r) = (\6{x) -6{x + r)\ p ) ~ r c ". (2) 

The set of scaling exponents, ( p , fully characterizes intermittency. In particular, deviations from 
the dimensional (linear) behavior C,2 P = PC2 are evidence of a nontrivial scalar transfer among 
scales, in analogy to the energy intermittent cascade in turbulent flows (see Ref. @ for a recent 
review on this subject). 

Another important issue is connected to the questioned universal character of scalar fluctuations 
in the inertial range. Universality should be the consequence of local, or quasi-local, interactions 
among scales and it should be reflected in a strong robustness of the scaling exponents, ( p , against 
variation of forcing and/or dissipation mechanisms. This idea dates back to the phenomenological 
work of Obukhov and Corrsin 

Statistics of the passive scalar must be strongly related to the properties of the advecting 
velocity fields; therefore, the only realistic case would be to study the "full problem" given by ([!]) 
together with the corresponding Navier-Stokes equations, in the limit of high Reynolds number, 
which describes the time evolution of a fully developed turbulent velocity field. 

Due to the lack of knowledge about the statistical properties of a fully developed turbulent flow 
the "full problem" cannot yet be attacked. Nevertheless, equation (D) with a prescribed velocity 
field of suitable statistical properties, could be a good "training ground" for understanding general 
mechanisms (if any) underlying the presence of anomalous scaling in fluid-dynamics. 

Recently, much attention has been paid in this direction by studying the intermittency proper- 
ties of a passive scalar advected by a velocity field which is Gaussian, white-in-time and which has 
self-similar spatial correlations. In Refs. j|, [| 0. [], [T^, |TU| the two-point velocity field correlation 
function was set to: (v i(x, t)v j(x', t')) = 5(t — t')Dij(x — x'), with D(x) = -D(O) — D(x). Here, D 
is the ci-dimensional velocity-field structure function : 

&ij{x) = D \x\* (d — I + £)5ij — £,XiXj\x\~ 2 , (3) 

where the scaling exponent £ of the second order velocity structure function, £, with < £ < 2, 
is a free parameter. Higher order velocity-field correlation functions are fixed by the Gaussian 
assumption. 

Such a choice is far from being realistic because of the assumed fast-time-decorelation and 
Gaussianity. Therefore, no quantitative agreement with the intermittent properties of the "full 
problem" must be expected. Nevertheless, in Refs. [|, || |7|, || [10], [Ll| many interesting analytical 



and phenomenological results for this toy-model have been obtained, giving for the first time the 



impression that some new insight into the "intermittency problem" has been gained (see Ref. [[13 
for related results in a passive vector case). 

The big advantage of this toy- model is that, due to the delta-correlation in time, all equations 
for equal-time n-points correlations are formally closed, that is, the n-th-order correlation depends 



only on lower order correlations. This fact was stressed by Kraichnan j|, ||, who gave for the first 
time the closed expression for ^(r). In Ref. Q a theory for all structure functions was proposed 
and an explicit formula for ( p given. The main physical outcome was that all structure functions of 
order greater than two have intermittency corrections and that intermittency should be connected 
to some non trivial matching between advective and diffusive properties of the model. Indeed, 
for structure functions of order greater than two it is no more possible to clearly separate the 
contributions from diffusion and advection terms in the generation of anomalous scaling. 

In Refs. [||, |8|, r0| it has been shown that intermittency of the scalar structure functions is con- 



nected to the properties of the null-space of the linear operators describing the inertial-scale evo- 
lution of multipoint passive scalar moments. Moreover, perturbative expressions of intermittency 
corrections, as a function of the parameter £ and as a function of the inverse of the dimensionality 



of the system were given [0, RL [L0j]. In both cases careful analysis of matching conditions at infrared 



(IR) and ultraviolet (UV) scales must be taken into account. Universality in the scaling exponents 
is still preserved, while constants in front of the power laws acquire dependencies on the forcing 
and on the passive scalar dissipation. How far all these results can be freely extended to the "full 
problem" is difficult to say. Certainly, the scaling exponents have a strong quantitative dependence 
on the statistical properties of the prescribed velocity field. Influence of a short-but-non-zero cor- 
relation time was studied in Ref. and found important, while from a phenomenological point of 
view it is difficult to say if the perturbation expansion of Ref. @ can be extended up to the value 
£ = 2/3 which mimics the spatial energy spectrum of true turbulence. (Let us notice that, even for 
£ = 2/3, the assumed delta-correlation in time for the velocity makes the classical Obukhov-Corrsin 
theory inapplicable to the passive scalar two-point correlation function ||.) 

An important drawback of the Gaussian delta-correlated toy-model is that numerical simula- 
tions seem even more difficult than a direct numerical simulation of the "full problem". This is 
because (HD is a stochastic differential equations (SDE) the solution of which requires huge com- 
puter resources, because of the massive use of random number generators and also because of 
the paramount difficulties in the implementation of higher than first order discretization schemes 



22 1 . The state of the art of SDE numerical simulations is hardly compatible with the high spatial 



resolution needed to study scaling properties of the solution of (H) with a white-in-time velocity 
field (see however Ref. for a first attempt at overcoming these difficulties, albeit at the prize of 
introducing some anisotropy into the velocity field). 

In this spirit, we have tried to work with the simplest toy-model which has some connection 
to the physics of the "full problem" , while at the same time aiming for something more tractable 
analytically and/or numerically. We have thus investigated the intermittency properties of a shell- 
model for a passive scalar advected by a prescribed stochastic velocity field. Shell models (see 
Ref. [14]] for a pedagogical introduction) have been already successful in helping to understand 



many issues connected to fully developed turbulence [15], [L6L [L7|, [18], 0, 



The problem of defining a shell model for the advection of a passive scalar by a deterministic 



and chaotic velocity field has been already investigated in Ref. |21[ . Here, we are going to present 
our numerical and analytical results obtained in a weakly modified version of the shell model 
defined in Ref. [F2TJ, but using a delta-correlated surrogate for the advecting velocity field. The 



main drawback of our model is that, following shell- model philosophy, it retains only nearest and 
next-nearest neighbor interactions in the shell index. Having only semi-local interaction in the 
Fourier space could be a problem for describing the physics of passive fluctuations for those types 
of velocity- field-correlations which induce strong non-local effects. For example, simple dimensional 
considerations tell us that the limit of very small £ should be dominated by strong non-local effects 
(which leads to an ultraviolet divergence of the eddy-diffusivity in the limit £ ^ [[J ) . Therefore 
we do not expect that our shell model could properly mimic the physics of the PDE ([!]) in the 
range of parameters choice < ( < 1. Yet, we are quite confident that our model captures 
very well the passive scalar physics for £ = 0(1), that is when dominant interactions become 
more local in Fourier space, as is the case for true turbulence. The major advantage of using a 
shell-model is that, now, reliable numerical simulations become feasible, allowing nonperturbative 
determinations of scaling exponents up to high orders (in this paper we present results up to order 
16). Furthermore, having reliable numerical simulations should open the possibility of relaxing 
some of the constraints put on the statistical properties of the velocity field, thereby allowing 
investigation of increasingly realistic problems. 

We found that our shell model shares many properties with the original SDE (fj). Among 
them, the most importants are: (i) the second order structure function has normal scaling, (ii) all 
structure functions of order larger than two have anomalous corrections, (iii) there is a remarkable 
coupling between both UV and IR scales with inertial terms in the closed equation satisfied by 
fourth-order shell correlations. This result suggests that the anomalous scaling is due to nontrivial 
contributions of integral and diffusive scales on the inertial-range properties. Moreover, we find 
that anomalous behavior tends to vanish when approaching the laminar regime, £ = 2, for the 
advecting velocity field. 

Even though we have not obtained an analytical expression for the anomalous scaling expo- 
nents, we are confident that many useful insights can be obtained by further investigation of the 
present model, or even simpler one. For example, it should be possible to define some exactly 
solvable models, where the shell-velocity correlations is chosen such as to give exactly solvable 
(but nontrivial) linear operators describing the structure functions inertial properties 

The paper is organized as follows : in Section § we introduce the random passive scalar shell 
model and we discuss some problems connected with its stochastic differential equation (SDE) 
formulation. In Section ^ we explicitly derive the closed deterministic equations for second and 
fourth-order moments, which are numerically solved in Section |] to obtain the scaling exponents 
of second and fourth-order structure functions. Monte Carlo simulations of the SDE are presented 
in Section |5|. Concluding remarks and suggestions for further work are given in Section [|. 

2 The model and its stochastic differential equation for- 
mulation 



We recall the main features of the passive scalar shell model of Ref. . The model is defined in 
terms of a shell-discretization of the Fourier space in a set of wavenumber defined on a geometric 



progression k n = k 2 n . Passive increments at scale r n = k" 1 are described by a complex variable 
6 n (t). The time evolution is obtained according to the following criteria: (i) the linear term is 
a purely diffusive term given by —Kk^9 n ; (ii) the advection term is a combination of the form 
k n 9 n ru n "] (iii) interacting shells are restricted to nearest and next-nearest neighbors of n; (iv) 
in the absence of forcing and damping the model conserves the volume in the phase-space and 
the passive-energy E = J2 n \@n\ 2 - Properties (i), (ii) and and (iv) are valid also for the original 
equation ([!]) in the Fourier space, while property (iii) is an assumption of locality of interactions 
among modes, which is rather well founded as long as the power law spectra of the passive scalar 
E(k) ~ k~ a has 1 < a < 3. 



Our passive scalar shell model, inspired by Jensen, Paladin and Vulpiani's model |[21|| , is defined 
by the following equations (m = 1, 2, . . .) 

[± + K k 2 m )9 m (t) = fkldWdW-CiWCiW) 
+ MCi(t)C 2 (t) +C 2 W«m-iW) 

+ c m (9* m+2 (t)u m+1 (t) + 9* m+1 (t)u* m+2 (t))} + S ml f(t), (4) 
where the star denotes complex conjugation and 

7 _ nm _ ^ m 7 k m —\ ^ k m j^\ , . 

™m " i ®"in "q - ' 2 ' ^ m 2 ' 

with m_i = uq = 9-i = 9 = as boundary conditions. The forcing term 5i m f{t) acts only 
on the first shell. Note that (|5]) is one of possible choices for the parameters a m , b m and c m 
ensuring the conservation of passive scalar energy J2 m ®m9* m and phase-space volume. In numerical 
implementations, the model is truncated to a finite number of shells iV (here N = 19) with the 
additional boundary conditions un+i = un+2 = #tv+i = @n+2 = 0. 

Our model differs from that of Ref. |H[ by the absence of complex conjugation on two of the 
u factors in the rhs of (01). The reason for this change will become clear below. Furthermore, in 



Ref. [21] the passive scalar model (HJ) was coupled to the GOY shell model [15, 16, 17] for the 
nonlinear dynamics of the velocity variables. The GOY model displays multifractal behavior for 
the w n 's and an interesting intermittent behavior is also found for the passive scalar. 

Our goal, as in Kraichnan's work 0], is to use a nonintermittent velocity field and then to find 
if the passive scalar is nevertheless intermittent. For this, we assume that the velocity variables 
u m (t) and the forcing term f(t) are independent complex Gaussian and white-in-time, that is, 
delta-correlated. Furthermore, as in Ref. ||, we make a scaling assumption for the spectrum of 
the ii m (t)'s, namely 

(u m (t)u* m (t')) = 5(t-t')D m , (6) 

where 

D m — k m ^ (7) 



will be called the velocity spectrum. For f(t) we assume (f(t)f(t')) = 5(t — t'). 



As long as the velocity variables u m (t) have a finite correlation time and, hence, smooth sample 
paths, there is no particular difficulty in giving a meaning to the set (§) of random ordinary 
differential equations (ODE's). A well-known difficulty arises with ODE's having white-in-time 
coefficients : the mathematical meaning of the equation is ambiguous The physicist's and thus 



our view-point is to define the solution of such "stochastic differential equations" (SDE's) to be the 
limit, as the the correlation time tends to zero, of the solution of a random ODE with nonwhite 
(colored) coefficients. In the mathematical literature on stochastic differential equations this is 
called the Stratonovich version of stochastic calculus |]22|| . When numerically solving a SDE such 
as (|D one cannot use standard numerical methods which assume smooth coefficients. It would also 
be highly impractical to use a smooth approximation of the coefficients with a small correlation 
time, since this requires time steps much smaller than the correlation time. Fortunately, there 
is an alternative formulation of the SDE, the Ito version, which overcomes this difficulty. The 
solutions of SDE's are Markov diffusion process which can be characterized by their diffusion and 



drift coefficients (also called Fokker-Planck coefficients). The so-called Ito equation p2| encodes 
this information in a way which leads naturally to efficient numerical schemes. 

Let us briefly explain how this Ito equation is obtained. For this it is better to work with an 
abstract form of the starting equation (with the forcing term omitted) : 

AQ 1 t 

- = M e + -M 1 (-)6. (8) 

at e e 2 

Here, M is a deterministic operator and a random Gaussian operator with finite correlation 
time. The e-scaling is chosen in such away that, for e — > 0, the lSA 1 operator becomes delta- 
correlated in time. Just as ordinary white noise may be written as the derivative, in the sense of 
distributions, of the Brownian motion process, we can write 

1,, ,t. dW 

where W(t) is an operator-valued Brownian motion. The Ito SDE associated to (|8|) is 

d9 = (M + V)9 + dW9, (10) 

where 

/•CO 

V= / (Mi(s)Mi(0))ds (11) 

is called the drift operator. The presence of the drift term in ([H]) is easily understood : a Neumann 
expansion of 9(t + dt) — 9(t) in powers of Mi, starting from (g), gives a second-order term, which 
is 0(dt 2 ) for fixed e but becomes T>9dt when the limit e — > is taken. The drift term may also 
be obtained by using the Ito formula [^] or Gaussian integration by parts (see, e.g., Ref. 0, 



Section 4.1). The term dW# is called the diffusion term. Since Brownian motion has independent 
increments over nonoverlapping intervals, the operator dW, which is 0{\fd~i) and has zero mean, 
is independent of the past and, thus, of 9. It follows that the diffusion term has vanishing mean. 
Hence, the mean of 9 satisfies 

^ = (M + V)(9). (12) 



The Ito SDE associated to (f|) reads : 

d9 m {t) = ((d m b m+ i + 6 m C m _ 2 )-D m _l + ( — <2m c m-l + frm+2 c m) D m+ i 
+ a m-lb m D m -2 ~ fl m+l c m^m+2 — K ^m)^rn(t) dt 

+ i[om(0m+i(*)am-idW^_ 1 (t) - 0^ (t)a m+ idW^ +1 (t)) 

+ &m(C-i(0am-2rfW*_ 2 (i) + ^_ 2 (t)a m _idW m _i(0) 

+ c m (^ +2 (t)a m+1 ^ m+1 (t) + ^ +1 (t)a m+2 ^,; +2 (t))] 

+ <W<W)(t). (13) 

Here, a m = fc£/ 2 , D m = the W m (t)'s and W/(t) are independent identically distributed 

complex- valued Brownian motion functions, normalized in such a way that (| W m (t)| 2 ) = (| Wf(t) | 2 ) = 
t. It is noteworthy that the drift term in (|13D , the top two lines on the rhs, which involves the sum 
of the energies in the neighboring and next-neighboring shells, may be viewed as an eddy-diffusivity 
term. For Kraichnan's original equation ([I]), the Ito equation approach has a scale- independent 
eddy diffusivity ]25]. Here, it is proportional to k~^. The difference stems from the absence, in 



the shell model of direct interactions between widely separated scales. We also observe that, if we 
had used the original model of Ref. pl[ with complex conjugates on all the u factors, the drift 
operator would involve nondiagonal elements coupling different shells, a situation we avoided. 



3 Equations for the second and fourth-order moments 

In this and the following sections we are interested in the scaling behavior of the p-th-order 
structure functions : 

((W /2 ) oc k~<>>, (14) 

where ( p is called the scaling exponent of order p. If ( 2p = V (2 the structure functions are said to 
have a normal scaling. If p 2p = ( 2p — V C2 7^ the scaling of the structure function of order 2p is 
said to be anomalous. 

It is well known that from a linear stochastic differential equation with white-noise coefficients 
it is possible to obtain exact equations for moments of arbitrary order 0, 0. For example, from the 
abstract equation @ one derives the closed equation for the first-order moment (|i~2|) . Higher order 
quantities such as9®9,9®9®9®9, ... satisfy also linear stochastic differential equations with 
white-noise coefficients, from which closed equations can be obtained for (9 ® 9), {9 ® 9 £g> 9 <S> 9) , 
etc. . Such equations may also be obtained by use of Ito calculus J22| . 

In the shell-model context the moment equations become excessively cumbersome beyond order 
four. We have obtained the closed equations for 

E m = (9 m 9*J and ^ = (9^9^). (15) 

The general structure of these equations is as follows 

Ei = (-2 K k 2 m 5 l>m + A hm )E m + F t . (16) 

Plm = [~2K(kf + k 2 m )5l, n 5mj + Bi m;nj ]P nj + Gl m . (17) 



The equations are written in explicit form in the Appendix B. Here, we stress a few important 
properties. The forcing F m is restricted only to the first shell m = 1. The matrix A^ m is symmetric 
band diagonal with bandwidth five and has the following scaling law : 

A l+S , m+S = k 2 s ^A l>7n . (18) 

Similarly, 

^l-\-s;m+s,n+s;j+s ^ s ^Bi mn j (19) 

Straightforward scaling arguments indicate that ( |IBJ ) and ( |TTD may possess steady-state solu- 
tions with normal scaling. Such solutions have 

E m oc k 2 ~t = k-& (20) 

and 

Pmrn fc m ^ 4 (21) 

with ^4 = 2(2. In order to find what kind of scaling actually holds, we now resort to numerical 
solutions of the moment equations (|T6| ) and fllTD- Since these are deterministic equations they 
can be solved with high accuracy at relatively low cost. This is not the case of the Monte Carlo 
strategy of Section 5, which allows us however, to tackle structure functions of high orders. 

4 Calculation of second and fourth-order structure func- 
tions from the moment equations 



For the numerical solution of (jig ) and (|17D we used a 19-shell truncation with values of k chosen in 
such a way that the diffusive cutoff is well within the available range of shells. The solutions were 
obtained by time-marching until the steady state is reached. Such a steady state is necessarily 
stable. The scaling parameter £ was varied in the range 0.2 < £ < 2. 

We found that the second order structure function E m displays always normal scaling with 

C 2 = 2-£ (22) 

the same value as obtained by Kraichnan || . To investigate the nature of the scaling of the fourth- 
order structure function P mm , we plotted it against the second-order structure function, following 
the, now standard, Extended Self Similarity (ESS) procedure (see, Ref. p3}] ). This gave us the 
anomalous part p 4 = £ 4 — 2( 2 of the scaling exponent, plotted in Fig. 1 against £. It is clearly seen 
that the scaling is anomalous (p 4 > 0), the anomaly is a decreasing function of £ and disappears 
as £ — > 2, the "laminar" limit, as in Ref. 0]. 

For small values of £, interactions become more and more nonlocal and a higher number of 
shells would be needed to have a sufficiently large inertial range. (The number of shells in the 
"intermediate inertial-dissipation range", in which the local Reynolds number is 0(1), grows like 



l/£.) Even though our shell model has only local interactions, purely kinematic effects introduce 
important long-range diffusive corrections in the limit £ — > 0. This is why we restricted our 
calculations of for £ ranging from 0.2 to 2.0. 

The fact that normal scaling is obtained for the second-order structure function is not very 
surprising, this being exactly the same situation as for Kraichnan's original problem. It is easily 
checked that the operator A appearing in ( |T6| ) has equipartition solutions (E m independent of 
m) in its null-space, both for the full operator and its iV-shell truncations. Such equipartition 
solutions have no associated passive scalar energy flux and cannot bring about anomalous scaling. 

Let us now consider the anomalous scaling for the fourth-order structure function which are 
connected to the properties of the operator B = Bi m ^ n j. The operator B has again equipartition so- 
lutions in its null-space, which cannot cause anomalous scaling. We have checked numerically that 
finite-shell truncations of B have no other eigenvectors in their null-space. This, only superficially 
contradicts the interpretation of anomalous scaling as arising from the zero modes of the inertial 
operator describing the evolution of moments || |10| . In our case there are no anomalous zero- 
modes of the operator B but still we see numerically a very clean intermittent behavior. Where 
does the observed anomaly come from? One can imagine two interpretations. The first one rests 
on the physical observation that the operator B is naturally long-range : it mixes inertial scales 
with cut-offs at both the UV and IR ends. Therefore, its inverse involves a nontrivial mixture of 
contributions from very different scales. This mixture could be the cause of anomalous scaling, 
thereby defeating naive (local) dimensional analysis. The second interpretation is to imagine that 
quasi-zero modes (which would become true zero-modes in the limit of an infinite inertial range) 
are already dominating the inverse of B. The two interpretations are not in contradiction. Indeed, 
our truncated system, being always influenced by UV and IR cut-offs, it naturally takes into ac- 
count boundary conditions and therefore never shows true zero-modes. A similar scenario takes 
place in Ref. 0, where zero modes in the infinite dimensional function space have to be matched 
with IR and UV physical cut-offs. 

Let us finally remark, that the analysis of the relevant eigenvectors in the eigenspace of the 
operator B is highly complicated due to the fact that we are looking at the cone of positive 
functions, which is not a linear space. 

5 Monte Carlo simulations for structure functions 

As we noted, the moment-equation strategy becomes inpractical for determining structure func- 
tions beyond the fourth order. We therefore resort to Monte Carlo simulations of the stochastic 
shell model in its Ito version ([13]). We used the "weak-order-one Euler" scheme, the details of 
which may be found in Appendix A. Roughly, this means interpreting the Ito equation (p~3|) as a 
time-difference equation. This scheme is of order one (in the time step At) for averaged quantities 
such as the structure functions. Averages are calculated as time averages, assuming ergodicity (we 
checked that changes in the seed of the random generator do not affect the results). Integrating 
over a large number of realizations is thus equivalent to integrating over many large-eddy turnover 
times. 



First we validated our Monte Carlo simulation by comparison with the results for the second 
and fourth-order structure functions obtained from the moment equations in the previous section. 
The comparison can be seen in Fig. 2. The agreement of the Monte Carlo and moment based 
calculations is comparable to the machine precision (single precision). We mention that, in order, 
to maximize the extent of the inertial range and avoid an inertial-diffusive range with algebraic fall- 



off of the passive scalar spectrum P3| , we assumed an exponential cutoff on the velocity spectrum. 

In all calculations 19 shells were used, the molecular diffusivity was varied between 2~ 10 and 
2 -8 and the time step between 2 -29 and 2~ 25 . In order to give an example of the quality of the 
scaling, we show in Fig. 3 the log-log plot of the 4 th , 6 th and 8 th structure function versus the 2 nd 
order structure function. This is done by ESS in order to improve the scaling. 

Fig. 4 shows the scaling exponents ( p determined by least square fits using ESS up to order 
p — 16 for £ = 2/3. The error bars are obtained from the least square fits. In all the calculations 
the inertial range included eight shells or more. When halving the number of samples in our 
statistics, we found that the values obtained remained well within the error bars. The graph of 
( p appears to be linear at values of p beyond 8, with a slope of about 0.29. We do not rule out 
that the asymptotic linear trend is an artifact due to insufficient statistics. We observe that the 
Cp's obtained with our Gaussian white-in-time velocity are less anomalous than those reported in 
Ref. PHI, where the velocity was already multifractal. 



6 Conclusions 

We proposed and studied a shell model for a randomly advected passive scalar. We suggest 
that anomalous scaling of structure functions of order greater than two is connected to nontrivial 
dependence of inertial-range properties on both integral and diffusive scales. This is shown, for 
example, in the closed equation for the fourth-order moments [Omd^diB*) , where inertial-range 
shell correlations depend explicitly on correlations between distant shells with m <C /. The main 
advantage of our model is that reliable numerical simulations become feasible. We numerically 
estimated anomalous exponents for structure functions up to order 16 and for various scaling 
exponents £ of the velocity field (0.2 < £ < 2). Comparing our numerical results with those 
obtained in a similar passive shell model |H] , but advected by a multifractal velocity field, we find 
large quantitative differences. This is due to the obvious fact that the ^-statistics are strongly 
correlated to the statistics of the advecting field. Nevertheless, the possibility of writing down 
closed equations for correlations of any order could help in the understanding of intermittency 
in more general cases. Unfortunately, we do not see how to implement in our model the kind 



of perturbation expansion done in Refs. |6|, |I0 |. The main difficulty is the lack of long-range 



interactions in our model which forbids a proper definition of an eddy diffusivity in the limit 
£ — > 0. Some long-range shell models should be introduced and studied if one wants to follow this 
path. 

Another interesting question is to find the simplest passive-scalar shell model which has inertial- 
range intermittency. For example, it is easy to extract from (^) an even simpler shell model, by 
retaining only terms with #m_ 2 (^) M m-i(£) an d Gm+2^) u m+iit) in the advecting part. This model has 



the same properties as listed in the Introduction and some preliminary numerical results indicate 
that it has similar intermittency corrections. 

Let us, finally, remark that, up to now, all the existing work on anomalies driven by stochastic 
velocity field was rooted in some kind of analysis performed in the physical space. In contrast, our 
shell model exists only in a kind of Fourier space. Understanding the anomalous scaling in our shell 
model at a phenomenological level could help in devising a phenomenology of intermittency. For 
example, ideas connected to the popular inertial-range cascade picture could be usefully revisited. 
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Appendices 

A The Euler scheme 



The weak-order-one Euler scheme [22] associated to the passive scalar shell model in its Ito form 
(0) reads : 

6>™ + = {1 + ((a m b m+ i + 6 m c m _ 2 )-D m -iAt + (— a m c m _i + At 

+ 0"m-lb m Dm-2^t — 0-rri+lCrnD m+ 2At — ft/c^ At) }6^ 

+ & m (Ci« m - 2 Aifr_ 2 + c- 2 « m -iA^i) 

+ C m (C + 2« m+ lA^ +1 + c + i«m +2 Awr +2 )] 

+ /5 m>1 AW;. (A.l) 



Here 



AWl = VAt V l, (A.2) 

where the rf^ are independent identically distributed complex random variables of the form a + ib 
where a and b are independent Bernoulli variables with values ±1/^/2. The Bernoulli variables 
are numerically generated by a linear feed-back shift register random number generator (see, e.g., 

The choice of Bernoulli variables rather than Gausssian variables is particularly convenient for 
numerical purposes. It ensures that averaged quantities such as moments are correct to first order 
in At, a choice consistent with the scheme. 



B Moment equations 

We give hereafter the detailed form of the equation for the second order moments E m = {9 m 9^) 
and the fourth-order moment P im = {6i9*9 m 9^} . 



-Ink't, + A m , m )E m + A. 



'-m,rn—2 J -'m—2 



E„ 



+ A 



m.m- 



-lE; 



m—1 



+ Am,m+lEm+l + -A 



rn,m+2E m +2 



+ F n 



:b.i: 



Pl,m — [1 — + kf) — A mm — Aii + ^-l,m2^ m +l, m + $1— 2,m^"^-m+2,m)]Pm,l 

A m m —2Pl,m—2 ^m,m—lPl,m—l ■^■m,m+xPl,m+l -^m,m+2Pl,m+2 
+ ^l,l-2Pm,l~2 + Al,l-lPm,l-l + Al,l+lPm,l+l + ^U+2-Pm,Z+2 

^/,m2 \A-m,m— 2P rn,m— 2 A mm —\P mrn —\ ■^■m,m+lPm,m+l ^m,m+2-Pm,m+2) 

+ (Si >m E m + SijEi + 5i im 8ijEi)Fi. 



(B.2) 



The following notation has been used. 





2' 


(B.3) 


F t 


= V 1^1, 


(B.4) 




= 2^I? m _i, 


(B.5) 


A m ,m—1 


= 2(a 2 m D m+1 + b 2 m D m _ 2 ), 


(B.6) 


A m ,m+1 


— 2(a^D m _! + c 2 m D m+2 ), 


(B.7) 




= 2c^D m+ i, 


(B.8) 


A 


^m+2,m ^m+l,m ^-m— l,m A m _ 2,m- 


(B.9) 
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